Pressure-induced Topological Phase Transitions in Rock-salt Chalcogenides 
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By means of a comprehensive theoretical investigation, we show that external pressure can in- 
duce topological phase transitions in IV- VI semiconducting chalcogenides with rock-salt structure. 
These materials satisfy mirror symmetries that are needed to sustain topologically protected surface 
states, at variance with time-reversal symmetry responsible for gap less edge states in Z2 topological 
insulators. The band inversions at high-symmetry points in the Brillouin zone that are related by 
mirror symmetry, are brought about by an "asymmetric" hybridization between cation and anion 
sp orbit als. By working out the microscopic conditions to be fulfilled in order to maximize this 
hybridization, we identify materials in the rock-salt chalcogenide class that are prone to undergo a 
topological phase transition induced by pressure and/or alloying. Our model analysis is fully com- 
firmed by complementary advanced first-principles calculations and ab initio-hased tight-binding 
simulations. 

PACS numbers: 73.20.At, 71.20.Nr, 71.70.Ej 
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Semiconducting chalcogenides CA^ in which the cation 
C is an element of the group IV (Ge, Sn, Pb) and the 
anion A is an element of the group VI (S, Se, Te), repre- 
sent an attractive class of materials, due to their unique 
structural and electronic properties. The most interest- 
ing compounds in this class — SnTe, GeTe, and the 
lead chalcogenides Pb^ — display rock-salt structure [1] . 
Both SnTe and GeTe have been long known for their 
ferroelectric properties in the low-temperature distorted 
structure [2 , whereas lead chalcogenides possess poten- 
tial relevance for thermoelectric and optoelectronic appli- 
cations [3^. Their outstanding properties have been char- 
acterized by a variety of experimental techniques [3l H] . 
In parallel, a large number of theoretical investigations — 
carried out with different methods ranging from ab initio- 
based tight-binding to density-functional theory calcula- 
tions — have addressed their peculiar electronic structure 
[SUTOJ. Interest in this material class has been renewed 
because of relativistic effects that may prove relevant for 
future spintronic devices. For example, a giant Rashba 
effect has been predicted in the ferroelectric phase of bulk 
GeTe [11]. Furthermore, SnTe in its undistorted struc- 
ture has been predicted to be a topological crystalline in- 
sulator (TCI) [12 ; its spin-polarized surface states have 
been confirmed by photoelectron spectroscopy later on 
[13]. These observations suggest to investigate the pos- 
sibility of topological phase transitions and their micro- 
scopic conditions in this class of narrow band gap semi- 
conductors. 

In the Z2 class of topological insulators, time-reversal 
symmetry ensures topologically protected edge states [Ml 
[T5] . An odd number of band inversions distinguishes a 
Z2 topological insulator from a conventional band insu- 
lator: the former shows an odd number of Dirac cones 
pinned at time-reversal-invariant momenta (TRIM), as 
found for example in Bii-^^Sba^ and Bi2Se3 [141 |T5]. I^ 



contrast, crystal symmetries play a central role in the 
class of TCIs [ini [17]. From a microscopic point of view, 
necessary conditions for the appearance of metallic sur- 
face states on insulating bulk structures comprise gen- 
eral symmetry requirements, strong spin-orbit coupling 
(SOC), and band inversions in the bulk electronic struc- 
ture. Concerning TCIs of the face-centered-cubic (fee, 
fm3m) space group, a mirror symmetry causes the ap- 
pearance of an even number of Dirac cones on (001) 
and (111) surfaces (which preserve the symmetry; see 
Fig. [ik and b). These cones are situated off the TRIMs 
F and X within the associated high-symmetry line of the 
surface Brillouin zone [12]. Since all five members of 
the cubic IV- VI class (i.e. GeTe, SnTe, PbS, PbSe and 
PbTe) belong to the fee space group, all of them and/or 
their alloys may be prone to topological phase transitions. 
Such phase transitions, characterized by a nonzero mirror 
Chern number \l2i HE] , have been proposed for PbTe and 
PbSe under pressure [l0l[T2]; moreover, they were experi- 
mentally observed in ternary alloys Fhi-x^UxA {A = Se, 
Te) as a function of doping fTQ'JS?. These findings are 
not simply explained by the strength of SOC: although 
SOC in Pb (Z = 82) is larger than in Sn (Z = 50), PbTe 
is nonetheless a conventional band insulator, as opposed 
to the TCI SnTe. This clearly calls for a deeper under- 
standing of the microscopic origin of the band inversions. 
In this Letter, by means of a semi-empirical analysis, 
we show that all members of the rock-salt chalcogenide 
class under consideration can be turned into TCIs under 
pressure, with the only exception of GeTe. The methods 
used comprise advanced first-principles simulations and 
ab initio-hdised tight-binding (TB) calculations. Consis- 
tent with sophisticated GW computations [T0|, we find 
that the fundamental gap of lead chalcogenides shrinks 
upon applying external pressure, closes at a critical pres- 
sure, and subsequently re-opens with an inverted band 



a) 

c 

A 



imm 






C) \sp.±l/2) 

^•«'>''^ \p,±3/2) \sp,±3/2) 




M 



|a±-?/2> 



|-yA±-^/2> 



FIG. 1: (Color online) Crystal structure of rock-salt chalco- 
genides CA (a) and the fee Brillouin zone (b), showing one 
mirror plane (green) containing the high-symmetry points 
r, L, and X as well as their projection onto the [001] sur- 
face. Schematic representation of the band structure at the L 
point in the topologically trivial phase (c) and upon pressure- 
induced enhancement of sp hybridization (d). Dashed lines 
separate occupied from unoccupied energy levels. 



character. This behavior is typical of a topological phase 
transition and exemplified for PbTe in Fig. [2] Here, we 
used a relativistic tight-binding scheme [6 for the semi- 
infinite system [2^ |24] . As long as the fee structure is 
preserved, these band inversions cause metallic surface 
states (Fig. (2)3). 

To benchmark the reliability of our predictions, we 
performed accurate density-functional calculations with 
hybrid functionals [25 , as implemented in VASP [26l 
[27] . These improve significantly with respect to the 
local-density approximation (LDA) or the (semi-local) 
generalized-gradient approximation, especially for nar- 
row band-gap semiconductors and lead chalcogenides [9] . 
For PbTe, we find the fundamental band gap at L to close 
at a volume ratio of ^/vo = 0.91 {Vq equilibrium volume; 
see also Fig. [2b) that is accompanied by band inversion 
(Fig.pl. This finding corroborates the topological phase 
transition deduced from tight-binding calculations (cfr. 
Fig. [2]). At variance with PbTe, hole pockets appear in 
the electronic structure of GeTe, concomitantly with the 
pressure-induced closure of the fundamental gap, thus 
triggering a semi-metallic state instead of a TCI. These 
essential differences require clarification in a microscopic 
picture. 

Assuming the undistorted rock-salt fee structure, the 
fundamental band gaps of CA compounds are located at 
eight equivalent L points in the Brillouin zone, at odds 
with many III-V and II- VI semiconductors in which the 
direct band gap is at F [8]. Many peculiar electronic 
properties of lead chalcogenides FhA were resolved by 
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FIG. 2: (Color online) Topological phase transition in PbTe 
upon pressure, calculated by a tight-binding method. For 
selected volume ratios ^/vo, the spectral density of the top- 
most layer of the (001) surface is represented in color scale (b, 
in states/eV). Topological trivial phases with mirror Chern 
number appear for the equilibrium volume (a) and close to 
the critical pressure (b). A topological nontrivial phase (c) 
shows a mirror Chern number of —2. Its Dirac surface state 
is depicted in a constant-energy cut around X at ^f — 0.05 eV 
(inset in c). A k path shown also in (c) is marked by arrowed 
lines. The k axis common to (a-c) shows 2/5 of the M-X and 
X-F lines around X (cfr. Fig.llp). 



noting that the Pb-65 band lies below the top of the 
valence band, thus forming valence bands with the p 
electrons of the A anion. The existence of an occupied 
cation-s band leads to strong level repulsion at the L 
point. Due to its symmetry, this explains the narrow 
band gap Eg and its unusual ordering within the series: 
Eg(PhS) > ^g(PbTe) > Eg(PhSe) [i[9]. LDA calcula- 
tions also predicted a negative pressure dependence of the 
band gap, which would suggest an inverted band struc- 
ture at L. However, GW calculations have shown that 
this inversion is a spurious outcome of the local approx- 
imation for the exchange-correlation potential, which is 
known to place the conduction band too low in energy 
(that is, too small a band gap) ^lOj . 



PbTe V/Vo = 0.91 
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FIG. 3: (Color online) Electronic structures of selected rock-salt chalcogenides at different pressures (indicated by volume ratios 
^/vo), with band characters highlighted by color scale (anions: blue; cations: red). Simulations are performed within DFT 
complemented by hybrid functionals. 



The central role of s electrons in determining the elec- 
tronic properties of CA compounds is deduced within 
the framework of linear combinations of atomic or- 
bitals ^8 . The only nonzero matrix elements at the L 
points in the nearest-neighbor approximation are those 
describing the hybridization between cation (anion) s 
and anion (cation) p states. Because of SOC, this hy- 
bridization close to the Fermi energy involves mainly 
a combination of atomic-like p states on C and A 
sites, with total angular momentum j = 1/2, namely 
I ±1/2) = (b.,TV2)±ib,,TV2)+|p.,±V2))/V^, and 
energies Cp^AiC) = ^p,AiC) ± 2A^(c); here, Cp^AiC) and 
^A(C) are the orbital energy and the spin-orbit coupling 
constant of the C {A) ion, respectively. In the atomic 
limit that corresponds to the topologically trivial phase, 
Ao = Cp^c — ^p,A > means occupied A and unoccu- 
pied C p shells (Fig. ^). The effect of the hybridiza- 
tion with the s states is to push both C and A p levels 
to higher energies; this energy shift is roughly propor- 
tional to the square of the effective sp hopping inter- 
action tsp and inversely proportional to the energy dif- 
ferences Ai = Cp^c — ^s,A and A2 = Cp^A — ^s.c- The 
energy gap at L can then be approximated as Eg ^ 



Ao + lOt^^ (^r^ ~ ^2^^) /3- Hence, for band inversion 
to occur, two conditions have to be satisfied: 1. the en- 
ergy separation between atomic-like p states (including 
SOC) must be sufficiently small, and 2. the sp hybridiza- 
tion should be strongly "asymmetric". This asymmetry 
is realized if Ai ^ A2, i.e. if the cation (anion) s level 
is energetically close to (far from) the anion (cation) p 
states. 

The band gap Eg is directly related to the k • p Hamil- 
tonian ai L^ H = mcFz [SI [12]; (7;^ = ±1 corresponds to 
the p character on cation C and anion A sites, respec- 
tively. Furthermore, a negative m = Eg/ 2 implies that 
conduction and valence bands at L would respectively 
derive from A and C ions. A reversal of m in the pres- 
ence of the mirror symmetries in the fee structure would 
imply a topological phase transition with an associated 
change of the mirror Chern number |12j. Since Aq = 
ep,c — Cp^^ — 2(Ac + A^), a large SOC is required to ful- 
fill the first condition. The second condition, on the other 
hand, could in principle be controlled by alloying ternary 
solid solutions CxC'i-xA (along the path pursued in 
[T9H22j ) or CAxA'i-x^ since all energy differences are ex- 
pected to change: within the virtual-crystal approxima- 



TABLE I: Estimated relevant quantities from Harrison's [28] 
and Lent's ( 6 , in brackets) parameterization. Spin-orbit cou- 
pling constants are taken from atomic values, as estimated in 
[29j . Lent's parameterization has been obtained through a 
fitting procedure including also d states. Energy differences 
A (in eV); lattice constants ao (in A) taken from [l]. 
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PbTe 


0.71 


10.08 


4.32 


6.462 




(L72) 


(12.70) 


(7.16) 




SnTe 


L35 


10.71 


4.75 


6.327 




(0.95) 


(13.13) 


(6.69) 




PbS 


3.15 


13.77 


1.89 


5.936 




(3.92) 


(16.25) 


(5.59) 




PbSe 


2.09 


13.29 


2.95 


6.124 




(3.29) 


(15.75) 


(5.18) 




GeTe 


1.19 


10.55 


6.63 


6.009 




(0.28) 


(11.92) 


(8.51) 





tion (VGA) one has A^ = xAP + (1 -x)Af ^^^^'\ Be- 
cause the band inversion is proportional to the strength 
of the sp hybridization, a straightforward way to induce 
a topological phase transition is to apply external pres- 
sure, as this directly affects the sp hopping tsp ex d~'^ 

[TomiEH]. 

Guided by these considerations, we performed an em- 
pirical screening within the family of rock-salt chalco- 
genides, using both Harrison's [28 and Lent's [6 param- 
eterizations (Table [l|. Even if figures are different, trends 
are consistent within both parameterizations. We can 
therefore loosely identify three subclasses: (i) The first 
subclass — that comprises PbTe and SnTe — is charac- 
terized by a relatively small Aq (mainly due to the large 
Axe) and similar differences A^^ — A^^. In this respect, 
the main reason why SnTe is a TGI but PbTe is not, 
could be the smaller equilibrium volume of the former 
with respect to the latter. In turn, a relatively small 
pressure could tune the topological transition in PbTe 
[12] and the doping-dependent topological transition in 
Pbi-a^Sua^Te could be ascribed to a (chemical) pressure 
effect, (ii) The second subclass comprises PbS and PbSe 
and shows Aq's approximatively twice as large as those 
of the first subclass, due to the smaller A^. However, 
A2 is much smaller than Ai; it is thus very likely that 
a reduction of the lattice constant and the associated in- 
crease of tsp could result in band inversion (due to the 
strongly asymmetric sp hybridization) and, hence, in a 
topological phase transition, (iii) Eventually GeTe, be- 
longing to the third class, even though displaying a small 
Ao, does not seem to fulfill condition 2. Thus, a reopen- 
ing of the gap after a pressure- induced closure is unlikely 
and a transition to a metallic state (rather than toward 
a TGI) is expected upon applying pressure (Fig. [3|. 

To verify the band inversions at the L points, we 
have performed ab initio hybrid-functional calculations 



[25][32. The band inversion is highlighted for PbTe, 
PbS, PbSe, and GeTe at selected volumes (Fig.|3|. For 
PbTe, the closure of the band gap shows up at a criti- 
cal volume of 0.921^; hence, the band inversion has al- 
ready occurred at a volume ratio ^/vo = 0.91 (pressure 
of about 4 GPa; our TB calculations give the same crit- 
ical volume. Fig. |2|. It shows a predominant Te anion 
character in the conduction band at the L point; the neg- 
ative gap increases upon further decreasing the volume 
(V'^o = 0.86 in Fig. [3]). Indeed, the atomic and orbital 
resolved character of the valence band maximum (con- 
duction band minimum) at L shows a predominant Pb-5 
and Te-p (Ph-p and Te-s) contribution at the equilib- 
rium volume, whereas an opposite character appears as 
the volume is decreased below the critical value [33]. The 
same trend with pressure holds for both PbS and PbSe, 
with critical volumes of about 0.941^ and 0.96^0 (tight- 
binding: 0.91Vb and 0.96Vb), respectively, corresponding 
to a pressure of about 4.2 GPa (6.3 GPa) for PbS and 
2.4 GPa for PbSe. For GeTe, on the other hand, the gap 
closes at a critical volume of 0.9^0, corresponding to a 
rather large pressure of 5 GPa. This finding is consistent 
with the smaller difference between Ai and A2 as com- 
pared to those of PbTe. Furthermore, since the j = 1/2 
and j = 3/2 manifolds in the conduction bands are close in 
energy due to the small SOG in Ge, a further increase of 
pressure would push them both below the anion p states, 
thus turning GeTe to a semimetal (Table [l| and Fig. |3|. 

To provide further support for the pressure-induced 
topological phase transitions, we performed relativist ic 
tight-binding calculations for bulk PbTe, PbSe, PbS, and 
GeTe. The computation of the mirror Ghern number is 
done in the spirit of the spin Ghern number [3^ . Since 
the Bloch states are eigenstates of both the Hamiltonian 
and the mirror operator [31 , we separate the Bloch states 
into two categories with mirror eigenvalues ±i and cal- 
culate the Berry curvature for both. The integral of the 
Berry curvature over the intersection of the mirror plane 
with the Brillouin zone yields Ghern numbers n±i, from 
which the mirror Ghern number Cm = (n+i — n_i)/2 is 
obtained [18 . While GeTe is always topologically triv- 
ial, we find a mirror Ghern number of —2 for band- 
inverted PbTe, PbSe, and PbS. In summary, the numer- 
ical calculations corroborate the microscopic picture de- 
rived above. 

In conclusion, we have shown that a strongly asym- 
metric hybridization between cation (anion) s and anion 
(cation) p states, together with a sizeable strength of 
spin-orbit coupling, is a necessary condition for band in- 
versions to occur at the L points and the related topolog- 
ical crystalline state to arise in rock-salt chalcogenides. 
By performing a thorough analysis of pressure effects in 
the entire family of fee chalcogenides, we verify the topo- 
logical nature of the transition, as shown by nonzero mir- 
ror Ghern numbers and gapless edge states at (001) sur- 
faces. Furthermore, we suggest that lead chalcogenides 



with topologically trivial properties could be turned into 
a topologically nontrivial state upon a combination of 
applied pressure and alloying through anion substitution. 
For instance, in the rbSe^^Tei-a^ alloy one expects an 
effect of chemical pressure (ao(PbSe) <C ao(PbTe)) and 
of the band shifts. The latter "asymmetrize" the hop- 
ping interactions, so that for a doping of x = 0.2 one 
would get Ao = leV, Ai = 10.72 eV, A2 = 4.04 eV, and 
ao = 6.39 A in the VGA. Also on the basis of the smaller 
critical pressure predicted for PbSe, a topological tran- 
sition is therefore likely to occur at reasonable values of 
applied pressure. In this context, we remark that alloy- 
ing two topologically trivial insulators is a yet unexplored 
path to engineer a TCI. This scenario may be relevant 
for the experimental search of conducting edge states, as 
the carrier concentration, largely determined by the pres- 
ence of cation/anion vacancies, can be easily controlled in 
lead chalcogenides during the crystal growth [3]. On the 
other hand, a high concentration of cation vacancies is 
frequently found in binary SnTe, resulting in an undesir- 
able p-type degenerate conducting behavior, a drawback 
that could be overcome by means of FhA anion alloying. 
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